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We use the density matrix renormalization group (DMRG) to perform a detailed study of the 
critical properties of the two dimensional Q state Potts model, including the magnetization and 
energy-density profiles, bulk and surface critical exponents and the Casimir amplitudes. We apply 
symmetry breaking boundary conditions toaixoo strip and diagonalize the corresponding transfer 
matrix for a series of moderately large systems (L < 64) by the DMRG method. The numerically 
very accurate finite lattice results are then extrapolated by efficient sequence extrapolation tech- 
niques. The critical density profiles and the Casimir amplitudes are found to follow precisely the 
conformal predictions for Q = 2 and 3. Similarly, the bulk and surface critical exponents of the 
models are in very good agreement with the conformal and exact values: their accuracy has reached 
or even exceeded the accuracy of other available numerical methods. For the Q = 4 model both the 
profiles and the critical exponents show strong logarithmic corrections, which are also studied. 
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I. INTRODUCTION 



The density matrix renormalization group (DMRG) 
method by WhiteO has considerably enhanced our abili- 
ties to diagonalize numerically the Hamiltonians of one- 
dimensional quantum systems. Among the systems 
treated so far by_the DMRG method we mention Heisen- 
berg spin chains!, Heiseaberg laddersEl and strongly cor- 
related electron systems! Recently, the first steps have 
been made towards the generalization of the method to 
other directions, such as two dimensional (2d) quantum 
systemsu, non-Hermitian HamiltoniansO, etc. 

According to numerical observations the DMRG 
method has a fast convergence, if the ground state and 
the excited states of the Hamiltonian are well separated, 
especially when the ground state wave, function can be 
approximated as a product of tensorsQ. This happens, 
among others, at the valence bond solid point of the. gen- 
eralized antiferromagnetic 5 = 1 Heisenberg chairoB. On 
the other hand the convergence of the method becomes 
slower at the critical region. The origin of this slowing 
down is two-fold. First, the spectrum of the density ma- 
trix at the critical point becomes more dense and there- 
fore more states have to be kept to guarantee a given 
accuracy. Second, the ground state of the Hamiltonian 
is nearly degenerate in the critical region and therefore 
the computational time needed to perform the Lanczos 
diagonalization increases considerably. 

In the first application of the DMRG to the critical 
rcgioixl the couplings of the Hamiltonian were subject of 
a renormalization procedure and the critical exponents, 
obtained from the fixed-point transformation, were less 



accurate than those obtained by the standard renormal- 
ization group. Later studies—concentrated on two di- 
mensional classical systemsEj~t3, have demonstrated that 
the thermodynamic properties of the critical systems can 
be accurately obtained by the DMRG method. In a 2d 
classical system it is natural to consider the row-to-row 
transfer matrix and calculate its leading eigenvalues by 
the DMRG method, as done first by NishindlO. More 
recently, the DMRG method and Baxter's corner trans- 
fer matrixes have been combined in an efficient iterative 
algorithmic. 

In this article we use the DMRG method to study the 
critical properties of the two dimensional Q state Potts 
model, in particular we focus on the behavior of local 
densities, as the magnetization and the energy density, 
and on the Casimir amplitudes for L x oo strips with 
symmetry breaking boundary conditions. These quan- 
tities have been derived from conformal invariance, but 
to our knowledge, never computed directly from a lattice 
calculation. We also calculate the complete set of surface 
and bulk critical exponents and compare our results with 
the exact (Q = 2) and conjectured values (Q = 3, 4) for 
these exponents. 

In the usual DMRG procedure the thermodynamic 
limit is approached by simply increasing the size of the 
system; typically one considers quantum chains of lengths 
of the order of 10 3 , which are generated by many DMRG 
iterations. Quantities as order parameters and correla- 
tion lengths are calculated from such large lattices. The 
drawback of this approach is that for such large systems 
the accuracy of the DMRG method is not as good as 
for chains of length of one order of magnitude smaller, 
especially at the critical point, for the reasons discussed 
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above. On the -other hand in the traditional finite size 
scaling analysiscJ, which is based on a numerically ex- 
act diagonalization of the Hamiltonian by the Lanczos 
method, the size of the finite lattices are of the order of 
10, which is often too small value for an accurate deter- 
mination of critical exponents from finite size scaling. 

In the present study we combine the standard finite 
size scaling analysis with the computational power of the. 
DMRG method. Using White's finite system algorithmt£l 
we calculate quantities with very high accuracy for strips 
of moderate widths (L < 64). The advantage with re- 
spect to the exact diagonalization of small lattices is that 
the finite size scaling is performed on a wider range of sys- 
tem sizes, corrections to scaling have weaker effects and 
one can determine critical exponents with a high accu- 
racy. Through the finite system DMRG algorithm one 
is also able to calculate accurately not only bulk, but at 
the same time also surface exponents, obtaining thus a 
deeper knowledge of the critical properties of the system. 

The structure of the paper is the following. In Section 
U we describe the DMRG method and its application to 
the Q state Potts model. Our results about the critical 
density profiles are presented in Section |l|, whereas the 
critical exponents and the Casimir amplitudes are inves- 
tigated in Section [V. Finally, in Section ^ we discuss 
our results. 



Since the model - except of the case of the Q = 2 
(Ising model) - is not solved outside of the critical point 
the critical exponents are not known exactly for general 
Q. Their vakiGS|, however, are conjectured by cpaformal 
field theorvc3~E2l and by approximate mappingsc2l. 

In a L x oo strip the local densities, such as the mag- 
netization 



m L (l) 



Q(S(s ltk )) - 1 

Q-i 



(4) 



and the singular part of the energy density 

e L {l) = (<5(s;,fc - Si+i,fc)) - e , (5) 

depend on the position of the layers I — 1, 2 . . . L. 

In the transfer matrix formalism the elements of the 
row-to-row transfer matrix Tj, are labeled by the possi- 
ble states of a row in the system: i — 1, 2, ... Q L . The 
leading eigenvalue A^ of the transfer matrix 



T L \v l ) = Al \v l ), 
jives the free energy per spin as 

/l = -— lnAi. 



(6) 



(7) 



II. THE DMRG METHOD AND ITS 
APPLICATION TO THE POTTS MODEL 

The Q state Potts modclEI - due to many exact and 
conjectured results - is an important testing ground for 
different approximate theories and numerical methods in 
statistical physics. We consider aLxoo lattice labeled 
by I = 1, 2 . . . L and — oo < k < oo with an Hamiltonian 
written as 

L 

PH = -K }^ {6 (si t k - s*+i,fc) - 5 (si >k - s(, fc+ i)} (1) 

k 1=1 

in terms of a Q component spin variable s^k — 0, 1, . . . , 
Q — 1 (/3 denotes the inverse temperature). In the follow- 
ing we take (3 = 1 and consider the ferromagnetic model 
(K > 0). In this case one can obtain the phase transition 
point of the system 



and the energy per bond at the transition point 



eo 



1 



(2) 



(3) 



from a duality transformation^. For general Q the model 
is integrable at the transition point and the. transition is 
of second (first) order for Q < 4 {Q > 4)EE In this pa- 
per we restrict ourselves to the second order transition 
regime. 



If the dominant eigenvector \vi) is normalized, the 
squares of its elements v\(i) are equal to the probabili- 
ties of finding a row in a state i. Average quantities as 
(<5(s/,fc)) can be calculated then: 



(8) 



where s;(i) denotes the value of the Z-th spin for a row 
in the i-th configuration. In this way one derives mag- 
netization and energy density profiles from the dominant 
eigenvector \vi)- 

The problem of the transfer matrix approach is that a 
numerically exact solution of the eigenvalue equation (|T^) 
is restricted to small lattices, since the dimensionality of 
the spin space grows rapidly with the strip width L. Fi- 
nite size results derived form transfer matrix calculations 
can be extrapolated only for a limited range of values of 
L. 

The DMRG provides a very efficient algorithm for the 
construction of effective transfer matrices Tj, of large 
L x oo strips, which are generated iteratively from a 
transfer matrix of a strip of a small width (Lq) which 
can be diagonalized exactly. In the initial step a row of 
the LqXoo strip is divided into two parts A and B, which 
are labeled by variables 1 < ja < Na and 1 < js < Nb , 
such that i = (j'a,Jb) and NaNb = Q L ° ■ Typically we 
split the system into two subsystems of the same size, i.e. 
= Nb, although this is not strictly necessary. 
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A B 

FIG. 1. Schematic view of the DMRG algorithm for the 
construction of effective transfer matrices Tl (the arrow de- 
notes the transfer direction). Starting from a small system 
(Lo — 4 in the example) the strip width is enlarged by two 
spins at each iteration. New spins are added at the center 
of the strip and rectangles denote blocks containing several 
spins, but labeled by m states only. 

In the construction of a larger strip only part of the 
states of the A and B subsystems are kept, namely the 
eigenstates of the density matrices of A and B corre- 
sponding to the largest eigenvalues. The elements of the 
density matrix for the subsystem A are defined as follows: 

N B 

p [A) (l A ,k A ) = ^2 VL ° ( 1 aJb) v Lo (k A ,j B ) (9) 

3B = 1 

where the sum is extended to all possible states of the 
subsystem B. The density matrix eigenvalues u> r , which 
we order according to wi > L02 > ■ ■ ■ > ^n a , are equal to 
the probabilities of finding the subsystem A in the corre- 
sponding eigenvectors |Q r ) when the whole system is in a 
state \vl )- In the calculation the spin space of the sub- 
system A is truncated and one keeps m eigenvectors |f2 r ) 
corresponding to the largest eigenvalues. The truncation 
error, defined as: 

rn 

e=l-^w r (10) 
r=l 

gives an estimate of the accuracy of the procedure. For 
rn = N A (all states kept) obviously e = 0; typically the 
density matrix eigenvalues u r decrease rapidly with the 
index r, so good accuracy can be reached with a moder- 
ate number of states kept. The truncated subsystem A 
is enlarged by adding an extra Potts spin; the old sub- 
system plus the spin corresponds to a new subsystem A, 
labeled now by mQ states. An analogous procedure is 
followed for the subsystem B, whose spin space is first 
truncated by the selection of the dominant m eigenval- 
ues of the density matrix for the B part p^ B \ and then 



enlarged by an extra spin. Combining the two new sub- 
systems A and B together one obtains a transfer matrix 
Tlq+2 for a (Lo + 2) x 00 strip. The procedure is repeated 
until the strip of the wanted width has been generated: 

Tl — * T La+2 — ► T La+4: T L . (11) 




A B 

FIG. 2. Schematic view of White's finite system method. 
Only one of the two subsystems grows at each iteration: sub- 
system A in the first two iterations and subsystem B in the 
third. The spins sa and sb sweep trough the strip width; 
from the transfer matrix one calculates, for instance, the 
average magnetization of the surface spin si. 

A schematic view of the first iterations is shown in 
Fig. @. In each DMRG step one solves the eigenvalue 
equation: 

f L \v L ) = ~X L \v L ) (12) 

defined in a Q 2 m 2 x Q 2 m 2 space. The quantity Jl = 
— Isx(Xl)/L approximate the exact free energy /j, of the 
L x 00 strip and from the dominant eigenvector \vl) one 
can calculate the magnetizations of the spins s^^.k an( I 
s L/2+i,fc at the center of the strip (see Fig. [I]) and the 
energy density of the bond connecting them. 

The itecations in (|ll|) describe White's infinite system 
algorithmic. Repeating the procedure many times one 
generates transfer matrices of very large systems: at each 
iterations new spins are added at the center of the strip 
and the boundaries are pushed farther away from each 
other (see Fig. |l|). Notice also that we implement open 
boundary conditions, as commonly done for DMRG cal- 
culations since it-has been found that the accuracy is the 
best in this casal3. 

For the calculation of the whole profile one must per- 
form more DMRG iterations keeping the system size 
fixed: 

f L -> Xf' -> jf ] -> . . . -> , (13) 
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as shown schematically in Fig. || (we have taken q sa 4L 
in our calculations). In each of the iterations of (|ll|) a 
spin is added only to one of the two subsystems, say to 
A, following the scheme described above. For the part B 
one takes smaller subsystems generated in the previous 
DMRG iterations. This procedure is continued until the 
subsystem B contains Lq/2 spins only, i.e. is known ex- 
actly. Then in each of the following steps a spin is added 
to B until the subsystem A contains Lq/2 only. In the 
present calculation this procedure has been repeated a 
couple of times. 

In the iterations (|l^) the two spins sa and sb of Fig. 
[I] sweep through the whole strip width and from the se- 
quence \vl)i • • - one calculates the magneti- 
zation and energy density profiles across the strip. At the 
same time the free energies calculated from the eigenval- 
ues of the transfer matrices Tj 



iteration 



L of p3| ) decrease at each 

< fi^ L> getting .closer and closer to the ex- 
act free energy /l (it is knownQ that DMRG free energies 
provide upper bounds for /x). 

The stops in ( |l3| ) describe the so-called finite system 
algorithmic, which is particularly indicated for the study 
of finite systems, since it provides higher accuracy than 
the infinite system method. The procedure is more time 
consuming, but has the advantage of providing informa- 
tion about the whole strip and it is necessary to calculate 
density profiles!!-! . 

Using the finite system method we have generated ef- 
fective transfer matrices for strips of widths L — 8, 16, 24 
. . . 64 at the system critical point. We have used m = 32 
for Q = 2, m = 60 for Q = 3 and m = 80 for Q = 4; in 
each cases the truncation error was e < 10 -10 . It is not 
possible to relate the value of e with the actual errors for 
the free energy and the density profiles. On the basis of 
our experience we expect an accuracy of 8 — 9 digits for 
the density profiles for L < 64 and a somewhat higher 
accuracy for the free energy. 



III. CRITICAL DENSITY PROFILES 



where (3, v and a are the magnetization, correlation 
length and specific heat critical exponents, respectively, 
while d denotes the dimensionality of the system. 
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In a critical system confined between two parallel 
plates, being a large but finite distance L apart, the local 
densities ($(?')) such as the magnetization and the energy 
density vary with the distance / as a smooth function of 
l/L. According to the scaling theory by Fisher and de 
GennesE^I 



(*(0)»6 



r 



F ab (l/L) 



(14) 



where ab denotes the boundary conditions (b.c.) applied 
at the two plates and x$ is the scaling dimension of the 
operator <i>. For the magnetization and the energy den- 
sity operators their value are connected with the more 
common critical exponents as: 

x m = f3/i' and x e = d — 1/v = (1 — cx)/v (15) 



FIG. 3. Scaled magnetization (a) and energy density (b) 
profiles for the Q = 3 Potts model with fixed and equal spins 
at the boundaries. The symbols corresponds to L x oo strips 
with L = 8 (circles), L — 16 (squares), L = 24 (diamonds), 
L = 32 (triangles up), L = 40 (triangles down), L = 48 (tri- 
angles left), L = 56 (triangles right) and L = 64 (pluses). 
The dotted lines are the conformal invariance predictions for 
the profiles. Also the values of the critical indices x m and x e 
used for the scaled profiles are indicated. 



The scaling function F a b in Eq. (14) has the following 
asymptotic behavior: 



F ab (l/L) = A 



1 + B nh ( - 



Z« 1 ' 



(16) 
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where B a b, the Fisher-de Gennes coefficients are con- 
nected to universal quantities. 

In two dimensions conformal invariance predicts the 
following form for ther-Oritical profiles with general con- 
formally invariant b.c.r^r 4 !: 



(*(0>ab 



L . 
— sin 

7T 



Gab 



(17) 



where the scaling function G ao {uj) depends on the uni- 
versality class of the model and on the type of b.c. ap- 
plied. Expanding Eq. (17) for (/L<1 one finds that the 
Fisher-de Gennes coefficients B a b in Eq. (16) are given 
by: 



B„ 



X$ _ 1 Ggb {I) 

6 2 G ab {l) 



(18) 
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FIG. 4. As in Fis. 1 but for Q = 4. 



In this section we present the critical profiles for the 
magnetization and the energy density calculated with the 
DMRG method and compare them with the predictions 
of conformal invariance (17). We have performed the cal- 
culations for (5 = 2,3 and 4 with all types of conformally 
invariant b.c. Here - for the sake of brevity - we present 
only the results for Q — 3 and 4. The Q = 2 (Ising) 
results, which can-be calculated more efficiently by the 
Pfaffian techniques, will be used in the next Section to 
calculate the critical exponents of the model. For all b.c. 
used we found that the Q = 2 magnetization and energy 
density profiles agree well with the conformal results of 
Eq. ©■ 

We start with the parallel spin b.c, i.e. we fix spins 
to the same state s = at the two boundaries and we 
indicate this choice by setting a p. b = in Eq. p?|) . 
According to conformal invarianceEH the scaling function 
in Eq. ([17]) both for the magnetization and the energy 
density is a constant: Goo = const. Notice that Eqs. (flj) 
and ( |17| ) imply that the scaled density profile L 1 * (4>(l)Jab 
depends on I only through the variable l/L. In Figs. || 
and ^ we show the scaled magnetization and energy den- 
sity profiles for the Q — 3 and Q = 4 models, respectively, 
together with the conformal results. As seen in the Fig- 
ures the DMRG data corresponding to different values of 
L collapse into single scaling curves, which are in good 
agreement with the conformal results. In the scaling plot 
we have used the values of the scaling dimensions con- 
jectured by conformal field theory and which are given 
in Table III. These scaling dimensions will be computed 



from a finite size scaling analysis of the DMRG data in 
the next Section. Notice that the magnetization profiles 
of Fig. ||(a) and Fig. ||(a) differ only slightly since the 
magnetic exponents in the Q = 3 and Q = 4 Potts model 
are very close to each other. A more visible difference 
between the Q = 3 and the Q — 4 case can be seen in 
the energy density profiles of Figs. (3|(b) and |](b). 

As next case we consider the fixed-free b.c, i.e. we fix 
the spins to the state s = only at one boundary of the 
system; we indicate these b.c. by setting a = 0, b = f 
in Eq. (|l7j). The scaling behavior of the critical densi- 
ties in a free, surface are governed by the surface scaling 
dimensionsES, which are generally different from the cor- 
responding bulk quantities. For the Q state Potts model 
the surface magnetization scaling dimensions 



(19) 



where [3 S is the surf ace magnetization critical exponent, 
are lis 



ad in Table III 



thcor 



as derived by conformal field 
We recall that the value of the surface energy 
scaling dimension is universally x s e = 2 for— two dimen- 
sional modelscj. The conformal prediction^ about the 
scaling function in the critical profiles in Eq. ( |l7| ) 
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for the magnetization and 



for the energy density {A and B are constants). 
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FIG. 5. Scaled magnetizations (a) and energy densities (b) 
for Q = 3 in strips with fixed free boundary conditions. The 
symbols are as in Fig. |^. 

Figures ^| and |^ show the scaled magnetization and en- 
ergy density profiles with fixed-free b.c. for Q — 3 and 
Q = 4, respectively. The data collapse into single scal- 
ing curves which are in good agreement with the confor- 
mal invariance predictions (|2^) and (^lj), except for the 
magnetization profiles of the Q = 4 state Potts model 
as shown in Fig. ^(a) where one finds strong deviation 
of the numerical data from the conformal invariance re- 
sults. The origin of these discrepancies is the presence 
of very strong, logarithmic corrections to scaling in the 



Q = 4 ruodefii As will be discussed in the next Section 
the leading 1 / In L corrections to finite size scaling are 
universal and they can be taken into account using effec- 
tive, size dependent scaling dimensions. We made use the 
same strategy here and for the largest finite system with 
L = 64 we compared the DMRG data with the conformal 
results, where in Eqs. ([l7]), ( p0| ) the corresponding effec- 
tive scaling dimensions were used. As seen in the inset 
of Fig. ^(a) we have obtained good agreement with the 
corrected formula. 
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FIG. 6. As in Fig. | for Q = 4. The inset of (a) shows the 
DMRG data for L — 64 (crosses) and the conformal profile 
of Eq. @, @ where we have used effective size dependent 
scaling dimensions x m (L) and xf n (L) given in ( pSl| ) and ( pi3| ) 
with L — 64. 
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FIG. 7. Scaled magnetizations (a) and energy density pro- 
files (b) in strip with mixed boundary conditions for Q = 3. 
The symbols are as in Fig. ^. 

Finally, we consider the third type of conformally in- 
variant b.c., when the two boundary layers are fixed to 
different states, say a = and 6=1. The conformal 
predictions about the scaling functions in Eq. (|l7]) are 
in this case the most complicatedEl For the magnetiza- 
tion profile the scaling function Gq\ (ui) is expressed by a 
combination of hypergeometric functions, which reduces 
to polynomials for Q = 2, where Gq^(w) = Aui and for 
Q = 4, where G^(w) = 3(1 + Alo + uj 2 ). For the energy 
density profile the scaling function 



Gq X — C 



l-xj2 
i — sii 

5 - 4x £ 



nl 



involves the energy scaling dimension x e 
denote non-universal constants). 



(22) 



(A, B and C 



FIG. 8. As in Fig. ^ for Q = 4. The inset of (b) shows 
the DMRG data for L = 64 (pluses) and the profile obtained 
from conformal invariance where we have used an effective size 
dependent critical index x e (L) given in Eq. (p2), for L = 64. 



Figures ^ and |g show the scaled magnetization and 
energy density profiles for Q = 3 and Q — 4 for mixed 
boundary conditions. In the Q = 3 case the numerical 
data referring to strips of different widths collapse into 
single curves which are in very good agreement with the 
conformal invariance predictions. For the case Q = 4 the 
agreement with conformal invariance is worse, again due 
to the logarithmic corrections. The strongest corrections 
are found for the energy density profile Fig. |^(b); the 
data collapse is clearly worse than in the Q = 3 case 
for the energy density. As done for the fixed-free bound- 
ary conditions we consider finite size approximants to the 
conformal profile (122), by substituting the exponent x e 
appearing in Eq. ( |22| ) with x e (L) for L — 64. Also in this 
case the modified conformal profile with an effective size 
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dependent exponent agrees quite well with the numerical 
data as seen in the inset of Fig. §(b). 



IV. CRITICAL EXPONENTS AND CASIMIR 
AMPLITUDES 

In this Section we perform a quantitative analysis of 
the finite lattice data obtained by the DMRG method. 
First, we consider the various critical exponents, then we 
investigate the finite size dependence of the free-energy 
density and calculate the Casimir amplitudes with differ- 
ent boundary conditions. 



A. Critical exponents 

The critical exponents such and x s m defined 

in Eqs. (15) and ( |l9|) can be calculated by two differ- 
ent methods. In the first method, which is more general 
and not restricted to two dimensional systems, we use 
the scaling form of Eq. ( |l4| ) . The second way of deter- 
mination of the critical exponents is based on conformal 
invariance: we use explicitly the functional form of the 



.(£) = 



density profiles presented in Sec. [II 



We start with the scaling method and compare the val- 
ues of a density, for instance the magnetization tulQ) for 
different widths L, which are taken as L — 8, 16,. . . 64, 
and at constant values of z = l/L. The corresponding 
scaling dimension - x rn - is obtained as the limiting value 
of the finite lattice estimates: 



\n(m L /m L+8 ) 
ln(l + 8/L) 



(23) 



which - according to scaling theory - is independent of 
< z < 1 (in the previous formula mj = m^zL) and 
m L+8 = m L+8 [z(L + 8)}). 

Following traditional finite size scaling methodsEj we 
extrapolate the finite lattice approximants x m (L) to L — ► 
oo, using a powerful sequence extrapolation technique, 
such as the widely used BST-methooH. Here, we are 
not going to recapitulate the details of this extrapolation 
procedure just for illustration we present the table of ex- 
trapolants (Table |) for x m obtained at the middle point 
(z = 1/2) of the magnetization profiles of the Q = 3 Potts 
model with parallel-spin boundary conditions. From the 
original series, which is in the first column of Table | and 
calculated through Eq. fl23|), new extrapolants are gen- 
erated by the repeated use of the BST algorithm, which 
are shown in the successive columns. One generally ex- 
pects faster convergence within the data in higher order 
columns. As one can see in Table | there is a nice conver- 
gence up to the third and fourth column, which however 
is not increased in the further steps. Thus our estimate 
of the magnetic scaling dimension based on the BST al- 
gorithm is 



0.1334(1) (Q = 3), 



(24) 



whiciL.agrees very well with the conformal result x m = 
2/15EJE9, and it is at least as accurate as the other ex- 
isting numerical resultsE3. 



TABLE I. Table of the BST extrapolants for the bulk magnetization exponent for the Q = 3 Potts model calculated from 
finite size scaling at the center of the strip with fixed identical spins at the boundaries. (The parameter of the BST algorithm 
is e = 5/6.) 



0.108817 0.128823 0.133723 0.133518 0.133505 0.133511 0.133715 

0.116774 0.131197 0.133577 0.133465 0.133510 0.133801 

0.120575 0.132097 0.133514 0.133431 0.133603 

0.122859 0.132535 0.133475 0.133409 

0.124400 0.132781 0.133449 

0.125517 0.132933 

0.126367 



For the calculation of the scaling dimension x e we have 
used an expression similar to Eq. ( p3|) with the magneti- 
zation replaced by the energy density profile. The surface 
exponent x s m follows from the analysis of the magnetiza- 
tion profile for fixed-free boundary conditions; according 
to scaling theory the magnetization of the surface spin at 
criticality in a L x oo strip scales as: 

m s L ~ L~< (25) 



from which one can extract the surface exponent in an 
analogous way as done for x m (see Eq. (|23|)). 

Next, we are going to describe our second way of cal- 
culating the critical exponents, which is based on the 
conformal results about the profiles. As an example we 
consider the magnetization profile with fixed-free b.c. in 
Eqs. (0), ( p0| ) and take three points from the profile 
for a strip of given width L: uil(L/4), uil(L/2) and 
mi(3L/4). It is easy to see that the surface magnetiza- 
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tion exponent x s m is approximated as: 

x s (L) = In [m L (|) /m L (f )] 
ml ' ln(V2 + l) 

whereas the bulk exponent x rn is involved in the combi- 



nation: 



(26) 



In 



mi 



In 2 



(27) 



TABLE II. Table of the BST extrapolants for the surface magnetization exponent for the Q = 3 Potts model calculated from 
Eq. (The parameter of the BST algorithm is e = 5/6.) 



0.554998 


0.673224 


0.668052 


0.667014 


0.666865 


0.584425 


0.670860 


0.667623 


0.666936 


0.666799 


0.600929 


0.669696 


0.667388 


0.666880 


0.666762 


0.611585 


0.669012 


0.667238 


0.666839 




0.619076 


0.668565 


0.667136 






0.624650 


0.668252 








0.628972 











0.666938 
0.666574 



0.666666 



Using other profiles one can easily obtain relations for 
the bulk (magnetization or energy) exponents, as well. 
To illustrate the accuracy of the conformal method we 
show in Table [n] the BST extrapolants of the surface 
magnetization exponent, calculated from Eq. (^6|) for 
the Q = 3 model. Here, one can see a similar tendency 
in the extrapolations as in Table [j] for the finite size data. 
The obtained estimate of the surface exponent 



x s m = 0.6667(1) 



(0 = 3) 



(28) 



is again in excellent agreement with Cardy's conformal 
result x s m = 2/3Eil. We note that our estimates on the 
critical exponents for a given model, obtained by the scal- 
ing and the conformal methods are consistent and the er- 
rors of the two methods are also comparable. The same 
is true if other type of profiles or other points of the pro- 
files are considered. We used this freedom to obtain an 
objective criterion about the error of our estimates. 

Our results for the critical exponents calculated by the 
DMRG method are presented in Table OH together with 



the exact and conformal results. The estimates are most 
accurate for the Q = 2 Ising model, which has two rea- 
sons. First, the structure of confluent singularities are the 
simplest for this model, and second, the numerical accu- 
racy of the DMRG method is also the most accurate in 
this case. The estimates for the Q — 3 model are still very 
accurate, as a matter of fact the accuracy of the DMRG 
method in this case is comparable or even highct-than 
other numerical, methods (Monte Carlo simulatiorS, se- 
ries expansions, traditional .finite size scalingo, corner 
transfer matrix approaclJijti3,etc). 

The Q — 4 state model, for which the estimates in 
Table III are less accurate, needs special considerations. 
It is known from a RG analysis that the critical Hamil- 
tonian of the model involves a marginal scaling field 'J, 
which results in logarithmic corrections to scaling. Then 
the finite-size scaling behavior of the free-energy density 
of the model as a function of^-t^j bulk (t, h) and surface 
(h s ) scaling fields is given by 



f(t, h, h a , L) = L- d f(L 2 - x °Z 3 /% L 2 ~ Xm Z x ' 1& h, L^^Z^h^Z^ , 1) 



(29) 



with 



Z 



*(o) 



InL 



(30) 



Now the different physical quantities can be obtained 
through differentiation and then their singular behavior 
at the critical point can be studied by finite size scaling. 
It is easy to see from Eq. (29) that the critical exponents 
determined by this way have universal 1/ In L corrections, 
which does not depend on ^(0): 



x m (L) = - 



1 1 
16hiL 



0[(lnL) 



^)4 + ^ + 0[(lnXn 



x s m {L) = l- — + 0[{\nL)- 2 ] 



(31) 



(32) 



(33) 



Then, the appropriate strategy, to analyze the finite 
size estimates for the exponentscll, such as for the bulk 
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magnetization scaling dimension x„ 
series of effective exponents x m (L) 
has then only (1/lnL) 2 corrections, 
type of analysis are presented in Table ttll which are 



, is to consider the 
- 1/(16 lni), which 
The results of this 



in satisfactory agreement with the conformal results and 
also give support to the validity of the RG analysis lead- 
ing to Eq. @. 



TABLE III. Comparison between exact (t) or conjectured (i") bulk and surface exponents for the Q < 4 Potts model and 
the numerical results obtained from finite size scaling extrapolations of DMRG data. 



Q 




x m (DMRG) 




x ms (DMRG) 




x e (DMRG) 


2 


1/8 T 


0.125000(2) 


1/2 T 


0.500000(5) 


It 


1.00000(1) 


3 


2/15* 


0.1334(1) 


2/3* 


0.6667(1) 


4/5* 


0.800(1) 


4 


1/8* 


0.120(5) 


1* 


1.02(2) 


1/2* 


0.49(2) 



TABLE IV. Comparison of Casimir amplitude predicted by conformal invariance (CI) and the numerical values obtained 
from extrapolations of the numerical data (DMRG) for different types of boundary conditions. 



Q=2 



Q=3 



Q=4 



A a0 (CI) 
,400 (DMRG) 
Ao f (CI) 
Aof (DMRG) 
Am (CI) 
Am (DMRG) 



-0.0654498 (-;§) 
-0.065447(5) 
0.13089969 (§) 
0.1309003(5) 
1.5053465 (^) 
1.505350(5) 



-0.10472 (-£) 

-0.104(1) 

0.287979 

0.2865(10) 

1.989675 (^) 

1.987(5) 



-0.1309 (-£) 
-0.127(5) 
0.6545 (§f) 
0.67(3) 
3.0107 (ff ) 
2.93(10) 



B. Casimir Amplitudes 

In the strip geometry with general conformally invari- 
ant boundary conditions the finite size dependence of the 
critical free energy-density is given as: 



fab(L) = f + 



fa + fb 
L 



L 2 



(34) 



where the bulk - /o - and the surface - f a + fb - contri- 
butions are non-universal. The second-order correction 
term A a b, known as Casimir amplitude, involves univer- 
sal quantities. As shown in RefsE-xtJ, the Casimir am- 
plitude and the Fisher-de Gennes parameter B a b in Eq. 
(pf) are related as 



Aab = - 



-B. 



ab ) 



(35) 



where c denotes the conformal anomaly number or cen- 
tral charge. We recall that c = 1/2 for Q = 2, c — 4/5 
for Q — 3 and c = 1 for Q = 4. From Eq. (p3i) combined 



with the results about the profiles in Section niM one ob- 
tains the values of the Casimir amplitudes predicted by 



conformal invariance. The simplest case is that of equal 
spins at the boundary for which the Casimir amplitude 
is just proportional to the conformal anomaly: 



Aoo = - Ta c. 



(36) 



In the numerical calculation we made use of Baxter's ex- 
act result about the bulk free energy-densitytS for the Q 
state Potts model which is given by the following formula: 



Jo 



lnQ 



dx sinh(7r 

— tanh(/ia;) — — - 

x smh(7rx) 



H)x 



(37) 



where cos/x = \/Q/2. 

Having subtracted /o from f a b(L) we extrapolated the 
surface contributions and then the Casimir ampltudcs 
using the BST algorithm. The results obtained for the 
Casimir amplitudes for different b.c. are presented in 
Table |y| together with the conformal results. As for 
the critical exponents the results for the Ising model 
(Q = 2) are the most accurate, but good agreement be- 
tween DMRG data and conformal results is found also 
for Q = 3. We stress that the numerical values for the 
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amplitudes are obtained from two successive BST extrap- 
olations. In the case Q = 4 also the Casimir amplitudes 
are affected by logarithmic corrections. From the analy- 
sis of the extrapolants for the conformal anomaly, which 
can be obtained numerically from Eq. (|36|), we found 
that such quantity has logarithmic corrections of type 
l/(lnL) 2 . There are no theoretical predictions available 
about the form of these corrections. In the case of peri- 
odic b.c. it is instead knownEa, from conformal results, 
that c has corrections of the type l/(lnL) 3 . 



V. DISCUSSION 

In this paper we have used the DMRG method to study 
the critical properties of the two dimensional Q state 
Potts model. Our study has demonstrated once again 
that, contrary to the widespread believe, the DMRG 
technique is able to describe accurately the critical region 
of a two dimensional classical (and consequently also of 
a one dimensional quantum) system. Indeed the numer- 
ical accuracy of the critical exponents calculated by our 
approach is comparable or even better than those of the 
other available numerical methods. 

As was noticed in early studies the DMRG technique 
faces into problems to describe the critical region, in the 
thermodynamic limit of the infinite system method. In 
this case, as was shown by Ostlund and Rommeru, the 
DMRG ground state is of a tensor product type, which 
can not describe algebraically decaying correlation func- 
tions, i.e. critical states. In the approach we used we 
have escaped this problem, since first, we applied the fi- 
nite system algorithm and we restricted ourselves to lat- 
tices of moderate witdh, and second, we applied symme- 
try breaking boundary conditions. In the finite systems 
studied the correlation length remains always finite (and 
the transfer matrix spectrum is gapped) , and the advan- 
tage of the boundary conditions is that one can calculate 
critical exponents from the non-vanishing value of the 
magnetization and energy density profiles. In the case of 
fixed-free boundary conditions from the same profile one 
can derive surface and bulk exponents. 

To analyze the critical point data we borrowed the 
method of the traditional finite size scaling and used effi- 
cient sequence extrapolation techniques, such as the BST 
method. As known from the theory of asymptotic series 
analysis, one can obtain accurate limiting value by the 
method if i) the terms of the series are numerically very 
accurate and ii) there are no strong confluent singulari- 
ties present. To satisfy the first requirement one needs a 
numerical accuracy of the data of at least 10~ 6 — 10~ 7 , 
which can be achieved by the DMRG method, even at 
the critical point, if moderately large finite systems are 
studied. These systems, as demonstrated by our present 
study, are still about one order of magnitude larger than 
those diagonalized by the Lanczos algorithm and used in 
the traditional finite size scaling analysis. 



In many physical problems the second requirement, i.e. 
the absence of (strong) confluent singularities is not satis- 
fied and these corrections to scaling contributions repre- 
sent the real limitations of the finite size scaling method. 
Since the strength of the confluent singularities decreases 
with the size of the systems one expects in many cases, 
like to our example with the Q = 3 and Q = 4 Potts mod- 
els, more accurate results about the critical exponents by 
the DMRG analysis compared with the traditional finite 
size scaling. 

In this paper we demonstrated another advantage of 
the DMRG method, namely at present it seems to be 
the only numerical method which could be used to study 
the density profiles in critical systems, especially in the 
parallel plate geometryM. Our numerical studies on the 
profiles, together with the evaluation of the critical expo- 
nents in Sec. have given strong numerical evidence 
that the conformal predictions about the critical density 
profiles of the Q state Potts model are exact in the scal- 
ing limit. Our results are satisfactory, even for the Q = 4 
state model, in which case the first order finite size loga- 
rithmic corrections were taken consistently into account. 

We mention that the analysis of this paper can be used 
and extended for other problems as well. For quantum 
Hamiltonians one expects similar or even better numer- 
ical accuracies, since the Hamiltonian is usually repre- 
sented by a more sparse matrix than the transfer ma- 
trix of the equivalent classical problem. By this method 
one can also study such systems in which the critical 
point is not known exactly by duality. Finally one could 
also study first order phase transitions by the present ap- 
proach, the results p£ this type of investigations will be 
presented elsewhereEl 
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